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Abstract 

We propose to model the uppermost inner core as an aggregate of randomly oriented 
anisotropic "patches" . A patch is defined as an assemblage of a possibly large num- 
ber of crystals with identically oriented crystallographic axes. This simple model 
accounts for the observed velocity isotropy of short period body waves, and of- 
fers a reasonable physical interpretation for the scatterers detected at the top of 
the inner core. From rigorous multiple scattering modeling of seismic wave prop- 
agation through the aggregate, we obtain strong constraints on both the size and 
the elastic constants of iron patches. In a first step, we study the phase velocity 
and scattering attenuation of aggregates composed of hexagonal and cubic crys- 
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tals, whose elastic constants have been pubhshed in the mineral physics literature. 
The predicted attenuations for P waves vary over two orders of magnitude. Our 
calculations demonstrate that scattering attenuation is extremely sensitive to the 
anisotropic properties of single crystals and offers an attractive way to discriminate 
among iron models with e.g. identical Voigt average speeds. When anisotropy of 
elastic constants is pronounced, we find that the S wavespeed in the aggregate can 
be as much as 15% lower than the Voigt average shear velocity of a single crystal. 
In a second step, we perform a systematic search for iron models compatible with 
measured seismic velocities and attenuations. An iron model is characterized by its 
symmetry (cubic or hexagonal), elastic constants, and patch size. Independent of 
the crystal symmetry, we infer a most likely size of patch of the order of 400 m. Re- 
cent hcc iron models from the literature are in very good agreement with the most 
probable elastic constants of cubic crystals found in our inversion. Our study (1) 
suggests that the presence of melt may not be required to explain the low shear 
wavespeeds in the inner core and (2) supports the recent experimental results on 
the stability of cubic iron in the inner core, at least in its upper part. 

Key words: multiple scattering, attenuation, dispersion, iron elastic properties, 
grain size. Earth's inner core 



1 Introduction 

The gross seismic features of the inner core are nowadays known with some 

confidence. The uppermost inner core, where crystallization of iron occurs, 

is a peculiar region with isotropic seismic velocities, strong attenuation and 
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where the presence of scatterers has been detected. Mineral physicists have 
concentrated a lot of efforts on experimental and theoretical investigations 
of elastic properties of iron at inner core condition. In spite of remarkable 
advances, fundamental questions pertaining to the symmetry class and the 
anisotropy parameters of iron are still actively debated and no clear consensus 
has emerged yet. In this study, we use the velocity and attenuation properties 
of seismic waves inferred from previous studies to propose a simple physical 
model of the uppermost inner core. Using a rigorous multiple scattering approach 
based on the Dyson equation for elastic waves in random media, we analyze 
critically various iron models from the mineral physics literature and give 
constraints on the possible stable iron phases in the inner core. In what 
follows, the necessary seismological and mineralogical data are reviewed and 
our approach is described. 



1.1 Seismic observations and mineral physics 



Traveltime analyses of seismic body waves have revealed that the Earth's inner 
core is anisotropic with about 1-3% velocity anisotropy and a fast direction 
of propagation parallel to Earth's rotation axis. Depth variat ions hav e been 



observed with an isotropic layer 



Su and Dziewonskil . 



1995 : 



Creager 



overlying deep anisotropy 



1999 



Garcia and Souriau 



Shearer . 



2000 



1994 



Wen and F. Niul . 



20021) and a change o 



( llshii and Dziewonskil . 



velocity anisotropy near the centre of the inner core 



2002 



Beghein and Trampert 



20031: 



Calvet et al. 



20061) 



The inner c ore is heterogeneous at differ e nt scales with hemispherical v elocit y 



variations (Niu and Wen. 



2001 



Garcia. 



continent scale ( jStroujkova and Cormier 



2002 



Cao and Romanowicz 



2004h . 



20041 ) and short wavelength features 
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(20-200 km) (IBreger et al. 



19991 ) . There are strong evidences for the concentration 
of very small scale heterogeneities at the top of the inner core. The PKIKP 
wave that bottoms in the uppermost inner core is strongly attenuated compared 
to outer core phases. In addition, its frequency spectrum is depleted in high 
frequencies. Various studies proposed a range of values for P-wave quality 
factor between 100 and 400 in the uppermost 200 km with possible lateral 



variations (e.g. 



2004 



Souriau and Roudil 



Yu and Wen 



1995 



Wen and F. Niu 



2002 



Cao and Romanowiczl . 



20061 ). V iscoelastic and scattering attenuation are cand idates 



to explain this observation (ILi and Cormier 



2002 



Cormier and Li 



2OO2I) but 



the latter mechanism seems more consistent with the strong coda of the 



reflected P wave generated at the inner core 



Poupinet and Kennett 




2004 


Levton and Kopei 




200' 


n. 1 



Koper et al 



boundarv (IVidale and Earle 



2004 



Krasnoshchekov et al. 



2000 



2005 



inner core are its isotropy and strong scattering properties. An important 
open question is to relate these macroscopic properties to the microscopic 
anisotropy of the iron alloy. 



Iron is thought to be the main constituent of Earth's inner core. However, 
the bulk properties and crystalline structure of Fe at such extreme physical 
conditions remain uncertain. The hexagonal-close-packed {hep) phase of iron 



has been proposed as a likely stable phase i n the inner core 



(e.g. 



-| r 



Stixrude and Cohenl . 119951 : Iyoo et al 



is no consensus on its anisotr opy (IMao et al 



1995 



j y seve ral authors 



Ma et al. 



1998 



2004h but there 



Steinle-Neumann et al. 



2001 



Antonangeli et al. 



the inner core is very large (j Williams et al. 



1994 



1993 



Ma et al 



. 1987: 


Boehler. 


1993: 


Saxena et al. 



20041). other structure s such as body-centered cubic ( 



Poirier and Shankland 



Saxena and Dubrovinsky 



19981 ) , double hexagonal closed packed (jSaxena et al. 
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19951 ) and orthorhombic (jAnderson and Duba 



1997 



Andrault et al. 



2000) 



cannot be excluded. The problem is made even more c omplex by the nature 



and quantity of light elements present in the inner core (jPoirier 



19M). Nickel 



or silicon may sta b ilise a body-centered-c u 



Lin et al. 



2002 



Vocadlo et al 



2003b 



jic (bcc) 



Vocadlo 



P 



lase at inner core conditions 



2007 



Dubrovinsky et al. 



20071 ). A long standing issue in mineral physics is the discrepancy between the 



Voigt average velocity of iron crystals and the observed se ismi c velocity. Even 



the m ost recent models such as proposed by 



Vocadlol fl2007h and 



Belonoshko et al. 



( 120071 ) ■ that incorporate the effect of temperature, have Voigt average shear 
velocity between 4.0 to 4.4 km.s~^. This is still significantly higher than the 
adopted 3.5 kni.s"^ velocity in seismological models. 



1.2 Approach of this study 



Two distinct classes of model have been proposed to explain the bulk seismic 
properties of the solid core: (1) Some authors have invoked the presence 
of elliptic liquid inclusions which could mimic the o bserved anisotropv and 



attenuati o n of s eismic waves through the inner core (jSingh et al. 



Bergman 



2000|); (2) 



(119971 ) proposed a solidification texturing model which results in 
depth dependent anisotropic properties in qualitative agreement with body 
wave traveltime measurements. In this work, we focus on the specific scattering 
properties of the uppermost inner core and its implication for the attenuation 
and velocity dispersion of seismic waves. Previous seismological studies of 
scattering have thus far only considered the effects of elastically isotropic 
heterogeneities with possibly anisotropic (anisomeric) spatial dis t ributions 



( jCormier et al. 



1998 



Vidale and Earle 



2000 



Cormier and Li 



2002 



Cormier 
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20071 ). These seismological models of heterogeneity are not clearly connected 
with the mineralogical and geodynamical models of the inner core. One may for 
instance ask: what is the physical meaning of an isotropic velocity perturbation 
in an intrinsically anisotropic material? 



We propose a simple model of solidification texturing of the uppermost inner 
core, consistent with seismological and mineralogical data. We limit our considerations 
to a few simple quantities that can be inferred from seismic observations: P and 
5" wave velocities and P-wave attenuation. We model th e superficial part of the 



solid core as an untextured aggregate of iron "patches" (IKrasnoshchekov et al. 



20051 ) as depicted in Figure 1. Each patch is characterized by the anisotropic 
properties of individual iron crystals and the orientation of crystallographic 
axes varies randomly from one patch to another. One patch is not necessarily 
made of a single crystal. One can for example imagine that it contains a large 
number of dendrites with strongly correlated orientation of crystallographic 
axes. Such a texture is shown in Figure 1 and is similar to w hat is observec . 



in laboratory experiments of crystalli zation of ice or hep iron (jWorster 



1997 



Bergman 



1997 



Bergman et al 



20031 ). The seismic properties of the aggregate 
depend on the patch size and the iron crystal properties (symmetry, anisotropy). 
As a consequence of our assumption of random patch orientation, the aggregate 
is macroscopically isotropic. But, because of crystal anisotropy, seismic wave 
velocities vary from one patch to an other which makes the aggregate microscopically 
inhomogeneous. A seismic pulse travelling through such a medium will be 
prone to amplitude attenuation and velocity dispersion. We show how multiple 
scattering theory can be used to relate the intrinsic elastic properties of iron 
crystals to the macroscopic seismic properties of the aggregate. An important 



conclusion of our theory is that the Voigt average velocity of a single crystal 
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is a poor approximation of the seismic velocity of the aggregate. 



2 Multiple scattering theory 



To describe statistically the aggregate, we introduce the following decomposition 
of the elastic tensor: 



ijkl 



Cijki - {Cijkii^)), 



(1) 

(2) 



where 



{6a,M) = 0. 



(3) 



Angular brackets represent the ensemble average over all possible orientation 
of a single crystal, also known as the Voigt average. The fluctuations 6Cijki{'x.) 
are considered to be single realizations drawn from an ensemble of random 
fields, having zero mean. In such an aggregate, the spatial variations of the 
elastic constants are only caused by different orientations of the patches. 
As usual the material heterogeneity is describ ed by the correlation of the 



fluctuations, which in our case is a tensor field (jStanke and Kino . 



19841 ): 



{SCijMix.)5Cai3y5i^')) 



^SlV|x-x'|), 



(4) 



where S denotes the eigth rank covaria nce tensor of the elastic moduli and 

i). The function 77(0;) 



1988 



depends on the symmetry of the crystal flHirsekornl . 
which equals 1 at x = and tends to as x — 00 is the spatial correlation 
function. It gives the probability that two points separated a distance x from 



one another are located within the same patch (jStanke 



19861 ). Equation 
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(4) reflects two assumptions regarding the statistics of the aggregate: flrst, 
there are no orientation correlations between different patches (i.e. there is no 
macroscopic anisotropy of the aggregate); second, the aggregate is statistically 
homogeneous. We have chosen a spatial correlation function of the form: 

ri{r) = e-""!". (5) 

The correlation function (5) implies that the patches (or grains) are convex 
and equi-axed, i.e., not elongated in some particular direction. Such a texture 
would produce anisotropic velocity and attenuation as shown by 



Margerin 



( I2OO6I ). In spite of the limitations mentioned above, an exponential function 



describes the variable shap es and linear d imensions of grains in a polycrystalline 

T). 



material reasonably well (iStanke 



19861 ). In this case, the effective average 



dimension of the patches is = 2 a. 

The inhomogeneity of the aggregate is related to the anisotropy of the patches. 
For example, if iron is weakly anisotropic, the elastic properties of the medium 
have only weak variations from patch to patch. The degree of inhomogeneity 



can b e expressed in terms of the effective elastic constants as (jStanke and Kino 



19841 ): 



- 4 (cy,)2 ■ 

To represent the fourth-order stiffness tensor, we have introduced in equation 
(6) the Voigt matrix Cu whose indices / and J vary from 1 to 6 with the 
following correspondence rule between tensor and matrix indices: (11) 
1, (22) ^ 2, (33) ^ 3, (23) 4, (13) ^ 5, (12) ^ 6. The degree of 
inhomogeneity is deflned with / = J = 3 for P waves, and with I = J = A for 
S waves. 
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To calculate the seismic respo nse of the aKgregate , we have used a formalism 



based on the Dyson equation (IRytov et al 



19891 ) for the ensemble averaged 



Green function. The theory takes into account all the physics of the problem: 
arbitrary anisotrop y of iron, mod e conversions between P and S waves, and 



multiple scattering (IWeaver 



19901) • This formalism yields effective phase velocities 
and spatial decay rates of coherent P and S waves propagating through the 
aggregate. At a given frequency u, the seismic coherent field (G) can be 



decomposed into longitudinal (P) and transverse (S) contributions (jWeaver 



1990): 



iG{u, p)) = (7^(p)pp + g^ip) (I - pp) 



(7) 



where g^{p) and g^{p) are defined as: 



(pv^n' + ^^'ip) 



uj^-{pVo^) +cjS{p) 



(9) 



In equation (7), I is the identity tensor, p is a unit vector in the direction of the 
wavevector p; Vff (Vq) and <7^{p) {a^{p)) denote the longitudinal (transverse) 
Voigt velocity and longitudinal (transverse) mass operator, respectively. The 
poles of the propagators g^{p) and g^{p) give the dispe rsion relation of P 



19951 ). To determine the 



and S waves in the aggregate, respectively (iSheng] . 
location of the poles, some approximation of the mass operators has to be 
made. In the First- Order- Smoothing App roximation, correct to second order 



in material heterogeneity. 



Weaverl fjl99Cll ) has shown that cr^{p) and a"^{p) 



are related to inner products of S with the unit tensor I and the second 
rank tensors pp. The detailed frequency dependence of the mass operators 
is governed by the correlation function rj. At sufficiently low frequency, the 
location of the poles can be given explicitly using the Born approximation 
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which consists in substituting p with ko = u/Vo in a{p) (iWeaverl . 119901 ). 
One obtains the following expression of the effective wave vector kg in the 
untextured aggregate: 



UJ 



a 



.Vr 



P/S 



1/2 



(10) 



When the Born approximation fails at high frequencies, one can still locate 
approximately the pole by calculating the density of states i n the r andom 



Shenej fll995h 



medium. This procedure is described in detail in the book of 
p. 85-86. In this work, Born approximation has been used to facilitate calculations 
and its validity has been verified a posteriori. 



The effective velocities V^^^ and attenuations a^^^ of the seismic waves are 
related to the real and imaginary parts of the mass operators, respectively: 



Ve 



P/S 



Vr 



P/S 



+ 



P/S 



P/S 



-s{a^'^{u/V,''")] 



:i2) 



Because the material is heterogeneous at the wavelength scale, the effective 
phase velocity in the aggregate is different from the Voigt velocity. The attenuations 
(y^P/s YQ^eci the amplitud e decay of t he co herent wave caused by scattering. 
In a recent review paper, iThompsonI (120021 ) shows the remarkable agreement 



between experimental measurements of attenuation in untextured metals and 
the theoretical prediction (12). To ease the comparison with seismic observations, 
we define P-wave and S'-wave attenuation quality factors as: 



Qp/s 



UJ 



(13) 
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3 Effective seismic properties of aggregates: study of 6 iron models 
from the literature 



3.1 Anisotropy and Voigt speed of single crystals 



Table 1 gives elastic properties of six iron crystals proposed for the inner 
core, obtained from either laboratory experiments or theoretical calculations. 
These crystals present very different anisotropic characteristics as revealed by 
Table 1. For example, the sign and amplitude of the anisotropic parameters 
for hexagonal iron defined in Appendix A vary significantly from one crystal 
to an other (Table 1). This is illustrated in Figure 2 where the P and S 
wave velocities of three hexagonal iron crystals are plotted as a function of 
the direction of propagation measured from the symmetry axis. In Figure 2a, 
the longitu d inal s ound velocity deduced from x-ray diffraction experiments 



flMao et al. 



19981 ) has a maximum at 45° from the symmetry axis. This 
result has not been supported by either first principle calculations or other 
experimental studies. The theoretical investigations themselves have led to 



contradictory 



conclusions: th e symmetry axis could be eithe r fast ( 



2000) or slow f Vocadlo 



20071 ) as shown in Figure 2. Recently, 



Laio et al 



Dubrovinsky et al. 



( 120071 ) have proposed that a body-centered-cubic (bcc) phase could be stable at 



inner core conditions. Re cent cubic iron models computed by f lBelonoshko et al. 



20071 ) and 



Vocadld (120071 ) at inner core temperature, present different anisotropic 
characteristics, described by the parameter u in Table 1 (Appendix A). All 
iron crystals in Table 1 have a higher S-wave Voigt velocity than that inferred 
from seismology -bet ween 3.5 km.s"^ and 3.67 km.s~^ in the akl35 model 



( iKennett et al. 



19951 ) -. Recent theoretical calculations of elastic properties of 
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iron at inner core conditions converge towards S-wave Voigt velocities between 
4.0 km.s~^ and 4.4 km.s~^. 



3.2 Seismic wave dispersion 



We explore the seismic properties of hep and bcc iron aggregates at the typical 
frequency of short period PKIKP waves (around 1 Hz). In what follows, 
an iron model will be simultaneously characterized by its single crystal and 
aggregate properties. We have calculated the P and S effective velocities for 
patch sizes ranging from 30 m to 100 km. As explained in section 2, the 
material heterogeneity induces a shift of the seismic wavespeeds from the Voigt 
velocities Vq. This effect is illustrated in Figure 3, where we have plotted the 
normalized variation of phase velocity 5V = (Vg — Vo)/Vo as a function of 
adimensional frequency ko a (where ko is the Voigt wavenumber and a is the 
correlation length) for hep aggregates (top) and bee agregates (bottom). Over 
the whole frequency range, the effective velocity of both P and S waves is 
reduced compared to the Voigt average. Typically, the shift is stronger by one 
order of magnitude for S wave than for P wave. For a given polarization the 
shift grossly increases with the degree of inhomogeneity (see equation 6), 
and varies over one order of magnitude at least, depending on the iron models. 
The effective velocity is also a sensitive function of the size of the patch. The 
strongest variations occur around k^^^a ~ 1 which corresponds to a strong 
coupling between P and S modes. The "jump" of the velocity dispersion curve 
around k^^^ a 1 is absent in scalar wave propagation an d is a remarkable 
featur e of the elastic propagation in metals as illustrated by 



Stanke and Kino 



( 119841 ). At very low frequency, our calculations verify the Hashin-Strikman 
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bounds as they should. Calculations with a gaussian two-point correlation 
function rj yield qualitatively very similar results. Since gaussian correlations 
are usually too smooth to represent natural materials, we have preferred the 
more realistic exponential function. 

Figure 3 illustrates the renormalization of the velocity in the heterogeneous 
material. Since the effective seismic velocity of an aggregate is lower than 
its Voigt average velocity, we can conclude that stable iron phases at inner 
core condition should have P-wave Voigt velocity close to the observed one, 
or slightly larger (around 2 %) if the degree of inhomogeneity is large. In the 
case of S waves, the difference between the Voigt and the se ismic velocity has 



to be positive, and can be large (up to 10 % for the model of 



Belonoshko et al. 



( 120071 )). Our calculations show that the Voigt velocity is not the appropriate 
parameter to discuss the compatibility of iron models obtained from mineral 
physics with seismic observations. 



3.3 Seismic wave attenuation 



We now examine the attenuation properties of six iron aggregates whose 
properties are summarized in Table 1. In Figure 4 we show the P-wave attenuation 
(l/Qp) as a function of adimensional frequency /c^a. Variations of l/Qp with 
respect to adimensional frequency depend on whether or a is varied. In our 
case, we work at constant central frequency u = 2tt and let a range from 
300 m to 90 km. A prominent feature is the strong frequency dependence of 
P-wave attenuation. Typically, the attenuation increases by about two orders 
of magnitude for P waves as k^a varies from 0.1 to 10. 
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Scattering attenuation is very sensitive to the anisotropic ch arac teristics of 



i ron c rystals. For example, the iron models of 



MaoetaL 



(119981) and 



Laio et al 



( I2OOOI ) have nearly equal Voigt velocities for both P and 5" waves (see Table 1). 
However, the elastic anisotropy of the crystals is quite different as illustrated 
in Figure 2. As a consequence, we find one order of magnitude difference for 
P-wave attenuation in the aggregates, at fixed patch size. At high frequency, we 
observe that longitudinal attenuation increases with the degree of inhomogeneity 
ep. All the attenuation curves show a smooth bump around k^a = 1, which 
corresponds to a frequency range where the coupling between P and S waves 
is maximum. Because attenuation and velocity are related by Kramers-Kronig 
relations, the "jump" of the velocity curve and the "bump" of the attenuation 
curve coincide. Considering these six iron models proposed for the inner core, 
we find as much as two orders of magnitude difference in attenuation for P 
waves at fixed patch size. Elastic wave attenuation is therefore an important 
quantity to characterize the properties of iron in the inner core. 



In Figure 4, the grey lines give some typical values for P-wave attenuation 



at 1 Hz at the t op o: 



by 



Yu and Wen 



the inner core, i.e. 300 < Qp < 600 as proposed 



(120061 ). Several iron models can explain the observed Qp 



by scattering attenuatio n only, but with different textures. For example, the 



Belonoshko et al. 



(I2OO7I ) model yields typical size of patches of abou t a few 



hundred meters, whereas the bcc iron model at 5500 K by 



Vocad 



Laio et al. 



ol fl2007h 



fl2000h . it 



implies patches around ten kilometers. In the case of 
seems difficult to explain P-wave attenuation in the uppermost inner core by 
scattering only. The study of these six iron models reveals a clear trade-off 
between elastic properties and typical size of patches (scattering objects). In 
what follows, this problem will be further examined in order to simultaneously 
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constrain crystal anisotropy and the texture of aggregates in the uppermost 
inner core. 



4 Constraints on elastic contants of iron in the inner core 



In this section, we consider the following inverse problem: find iron models 
that match the seismic measurements of P and S wave velocities, and P wave 
attenuation at 1 Hz in the uppermost inner core. An iron model is defined by 
the elastic tensor of a single crystal and the texture of the aggregate. Because 
the stable phase of iron at inner core conditions is still uncertain, we propose 
to examine two kinds of crystal symmetry: hexagonal and cubic. In our simple 
model, the texture is characterized by the correlation length a, equal to half 
the average patch size. We let a vary over more than two orders of magnitude: 
50 m < a < 10000 m. 

The observation of strong PKiKP coda suggests that scattering attenuation is 
the dominant mechanism for atte nuation in the inner cor e. T his interpretatio n 



mm and 



Cormier] (120071 ). 



is also supported by the study of I Cormier and Li 
A commonly accept ed value for Qp at the top of the inner core is about 300 



( Yu and Wen 



20061 ). We impose that an acceptable iron model should have a 
scattering Qp within the range 300-600. This choice implies that more than 
half the observed attenuation is explained by scattering and leaves some room 
for other attenuation mechanisms. Since S'-wave attenuation in the inner core 
around 1 Hz is still poorly known, we do not impose constraints on Qs- The 
range of acceptable seismic velocities corresponds to the values given by the 
akl35 model for the whole inner core. 
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4.1 Hexagonal iron 



To find acceptable iron models we perform a systematic grid search in the 
Cij space. Only positive definite elastic tensors have been retained. For each 
correlation length, we test 1.8x10^ physical iron models with anisotropy parameters 
(5,5,7) smaller than 50%. In Figure 5 (left), we represent the percentage of 
acceptable hexagonal models as a function of the correlation length. The curve 
shows a sharp maximum at a ~ 220 m, which yields an average size of the 
patches around 450 m. The correlation length of the aggregate cannot be 
smaller than 140 m but can be as large as 10000 m. It is to be noted that the 
most probable patch size is about 20 times smaller than the central wavelength 
of P waves. A careful analysis reveals that the bounds on the seismic velocities 
determine the overall number of acceptable models, whereas the bounds on 
the attenuation fix the position of the maximum. 

For three correlation lengths a =140 m, 220 m, 1000 m, we examine in Figure 
6 the distribution of acceptable models in the (e, 5) parameter space. We focus 
on the two parameters controlling the P-wave crystal anisotropy, in order to 
facilitate comparisons with results from mineral physics. Figure 6 is a 2-D 
histogram of the number of acceptable models. For each correlation length, 
there is a clear trade-off between e and 5 because only three independent 
observations serve to constrain 5 independent parameters. To avoid artefacts 
due to uneven sampling of the model space, we have taken care of normalizing 
the results by the total number of tested models per class. In fact the systematic 
exploration of the C/j space yields relatively homogeneous coverage of the (e, 
5) space. In Figure 7, we show the a-posteriori histograms of acceptable Voigt 
velocities for the three correlation lengths. 
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At small correlation length (a = 140 m), the aggregate is most likely composed 
of hexagonal iron crystals with a slow symmetry axis (20% < e < 50%), and 
possibly high Voigt velocities (11.6 km.s"^< V^f <12.1 km.s^^, 3.8 km.s~^< 
Vq <4.7 km.s~^). There is no clear constraint on 6. We remark that Voigt 
velocities of iron crystals can never be smaller than seismic velocities, but high 
Voigt velocities are not necessarily incompatible with seismic observations. 

At large correlation length (a > 500 m), the {e , 6) distribution has the 
shape of an ellipse centered around (0,0), with positive correlation between 
e and S in the range -10% < e < 10%, -20% < 6 < 25%. The crystal 
Voigt velocities must be close to seismic velocities. Careful analysis shows that 
there are slightly more hep crystals with a slow symmetry axis, than with a 
fast symmetry axis. At large correlation length, none of the 8 hexagonal iron 
models considered in this study has both anisotropic parameters and Voigt 
velocities compatible with seismic observations. 

For a ~ 220 m, we obtain a maximum of acceptable iron aggregates. Corresponding 
iron crystals have a complex distribution in e and 6. 6 is not constrained, and 
e varies between —20% and +50%. Interestingly, the elliptic domain which 
defined acceptable models at large correlation length is forbidden when a 
equals the most probable patch size. Like for small correlation lengths, we 
obtain a larger number of hexagonal iron crystals with a slow symmetry axis 
for P waves. The distributions of acceptable P-wave Voigt velocity (11.3 km.s^^< 
Vff <12.0 km.s~^) and S-wave Voigt velocity (3.5 km.s~^< Vq <5.5 km.s~^) 
are very broad, with an ave rage clearly higher than the seismic velocities. 



The iron crystal proposed by 



Steinle-Neumann et al. 



(120011 ) is the only model 



which presents both anisotropic parameters {e, 6) and Voigt velocities that are 
compatible with seismic observations. However, this iron model is questionable 
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because of t he too large value of the axial ratio of hep iron predicted by the 
calculations (IGannarelli et al.l . l2005l ) . 



4.2 Cubic iron 

We perform a grid search in the (Cn, C12, C44)/p parameter space, where p 
is the model iron density. The normalization by p facilitates the comparison 
between models with different densities. In any case, the scattering properties 
of an aggregate depend only on Cjj/p. 

For a given texture, we test 8.8x10^ cubic iron crystals with a homogeneous 
distribution in the range: 0.1 < Cu/p < 0.2, 0.06 < C^/p < 0.12, < 
p < 0.04, where the C/j and p units are GPa and kg.m~^, respectively. In 
Figure 5 (right), we represent the percentage of acceptable cubic models as a 
function of the correlation length. The curve shows a maximum at a ~ 200 m, 
which yields an average size of the patches around 400 m. Therefore the most 
probable size of patch is a robust feature, independent of the symmetry of 
the crystal. Like in the hexagonal case, the correlation length of the aggregate 
cannot be smaller than 150 m but can be as large as 10000 m. The percentage 
of accepted models is by one order of magnitude smaller in the cubic case than 
in the hexagonal case. 

For three correlation lengths (150 m, 200 m and 500 m), we examine, in Figure 
8, the distribution of elastic constants of acceptable cubic crystals. At large 
correlation length, we observe two distinct families shown in grey and black in 
Figure 8 (middle-bottom), whereas a single family exists at correlation length 
smaller than 200 m. The black and white circles show the recent results of 



18 



Vocadld (120071 ) for cubic iron at two different inner core tempe r ature s, and 



Belonoshko et al 



(120071 ) ■ This 



the black square is the recent bcc iron model by 
last model is completely compatible with seismic observations and corresponds 
to the "grey" family of cubic models. 

In Figure 9, we show the a posteriori distributions of Voigt velocities for three 
correlation lengths. The main features are very similar to the hexagonal case 
but the bounds are slightly sharper. For the most probable patch size, the 
average ^- wave Voig t veloc ity is around 4.4 km.s~^ , in ag reement with recent 



results by 



Vocadlol (120071 ) and 



Belonoshko et al. 



(120071 ). Our study shows 



that some cubic iron crystal with S-wave Voigt velocity around 4.0 km.s ^ 
can explain seismological d ata without invok ii ig the presence of melt, as for 



example the iron model by 



Belonoshko et al 



bcc iron model at 5500K by 



(120071 ). On the other hand the 



Vocadlol (120071 ) ■ although in reasonably good 



agreement with the results of our inversion, is not anisotropic enough to induce 
a strong renormalization of the seismic velocities in the aggregate (only 1 % 
as shown in Figure 3). For such an iron model, it is necessary to invoke the 
presence of some melt to explain the seismic velocities. 



5 Discussion and conclusion 



5.1 Patch size 



[1) From mineral physics and geodynamics. Estimates of the size of iron 



crystals in the . 
to about 5 mm. 



larth's inner core vary from the entire volume, some 1200 km. 



Stixrude and CohenI (119951 ) have proposed that a single crystal 



of iron with anisotropy around 4 % could explain seismic traveltimes. The 
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smaller estimates stem from geodynamical considerations. Near the melting 
temperature of iron, a grain size of about 5 mm is nece s sary t o make the 



viscosity smaller than 10 Pa.s, as proposed by 



Buffett 



dynamic recr ystallisati on, 



around 5 m. 



Bergman! (11997 



Yoshida et al. 



(Il997l ). Assuming 



fll996l ) obtained a typical crystal size 



19981 ) suggested that observed attenuation and 
velocity anisotropy of the inner core may result from radially elongated columnar 
grains and estimated the columnar grain width around 200 m. This last study 



is in good agreement with our most probable patch size. 



(2) From sei smology. With s c atteri ng and a fabric interpretation of seismic 



attenuation, ICormier and Lil (120021 ) have proposed average scale length of 
isotropic heterogeneities around 10 km, and P-wave velocity perturbations 
around 8.4% from the centre of the inner core to about 1000 km radius. 



Vidale and Earlel (l2000h 



From the energy envelopes of PKiKP coda waves 
have obtained a scale length of about 2 km, and elastic moduli fluctuations of 
1.2% in the uppermost 300 km of the inner core. These results may reflect a 
trade-off between velocity perturbations and scale length. However, the most 
serious criticism is the lack of physical connection between the hypothesized 
isotropic seismic heterogeneities and the crystalline nature of the inner core. 
Our most pr obable correlation leng th (a round 200m) is srnaller t han previous 



estimates by 



Cormier and Lil (120021) and 



to the lower limit proposed by 



Cormier et al. 



Vidale and Earld (120001 ) . but is close 



19981 ) (about 15 % P-wave 



perturbations and scale length in the range 0.5-2 km). 
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5.2 On the presence of melt 



The vast majority of proposed iron models for the inner core displays much 
higher Voigt average velocities than what is inferred from seismic traveltimes. 
This discrepancy has sometimes b een used as an argument to supp ort the 



presence of liquid in the inner core (jSingh et al. 



2000 



Vocadlo 



2007h . In our 



study, we have shown that high Voig t velocities, typica l 



is found in the most recent studies ( Belonoshko et al. 



y of the order o 



2007; 



Vocadlo 



what 



20071) 



do not necessarily imply the presence of melt. The reason is that the effective 
velocity in an aggregate can be as much as 15 % smaller (according to the 
anisotropic characteristics of iron crystal) than the average velocity in a single 
crystal. Two additional points have to be noted: (1) Iron crystals with Voigt 
velocities smaller than seismic wavespeeds have to be rejected, independent of 
our assumption for the value of Q (2) Iron models with Voigt velocities larger 
than 5.5km.s~^ can also be clearly rejected, unless they co ntain a huge a moun t 



of melt, typically more than 20 % according to Figure 2 in 



Singh et al. 



fl2000h 



5.3 Cubic vs Hexagonal 



Most theoretical and experimental studies show that hexagonal metals a t 



ambient temperature have a fast symmetry axis (lAntonangeli et al 



But t he effect of temperature on the elastic constants is still debated (ILaio et al 



2006h 



2000l ). We found that hep iron models with a slow symmetry axis are more 



probable in the inner core. If further experimental or theoretical investigations 
definitely show that hep iron at inner core conditions has a fast symmetry 
axis, we may conclude on the basis of our uppermost inner core model that 
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the stable iron phase is unhkely to be hexagonal. Recent cubic models from 
the literature are in very good agreement with the results of our inversion. 
Therefore, o ur study supports a stab le bcc iron phase in the inner core, as 



proposed by 



Dubrovinsky et al. 



fl2007h 



Although our study provides some indications on the possible elastic constants 
of iron, some large uncertainties remain. They are clearly caused by a lack of 
seismological constraints. In future studies, the modeling of the following three 
observations could significantly reduce the uncertainties on the properties 
of iron in the inner core: (1) energy envelope of PKiKP coda; (2) direct 
measurements of at 1 Hz; (3) velocity and attenuation anisotropy in the 
bulk of the inner core. As a large number of records of PKiKP coda are 
available, for which multiple scattering may be relevant, we plan to focus our 
efforts on this aspect in future works. The modelling of frequency dependent 
inner core attenuation combined with coda analysis could unravel the relative 
contributions of viscoelasticity and scattering attenuation. 
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A Anisotropic parameters for hexagonal and cubic crystals 



A.l Hexagonal symmetry 



In a transversely isotropic medium (hexagonal symmetry), the symmetry axis 
defines a particular direction which is usually chosen along the X Q-axis for 



conven ience. To describe anisotropic properties of hexagonal crystals. 



Mensch and Rasolofosaon 



e 



(I1997I ) have defined three anisotropic parameters: 

(Cii — C33) 

(C13 — C33 + 2C44) ( \ -[^ 

'-'33 

7 



2C44 

The parameter e represents P-wave anisotropy i.e. the difference between the 
phase velocities perpendicular and parallel to the symmetry axis. A negative 
value of e corresponds to a fast symmetry axis for P waves (Table 1). The 
parameter 7 represents the anisotropy of S" waves in a similar way. The parameter 
5 controls P wave propagation at intermediate angles from the symmetry axis. 
In order to fully prescribe the elastic properties of hexagonal crystals, two 
additional parameters have to be introduced. We have chosen the P and S 
Voigt average velocities because they are routinely compared with the seismic 
measurements. 
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A. 2 Cubic symmetry 



The stiffness tensor of a cubic crystal is described by tliree independant 
constants Cn, C12 and C44. The invariant anisotropy factor for cubic-class 
crystal is defined as: 

z/ = Cii - C12 - 2C44 (A.2) 

The elastic properties of cubic crystal are fully described by the parameter v 
and the P and S Voigt average velocities. 
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Hexagonal Iron Models 


P 


Cii 


Cl2 


C44 


C'ss 


Cl3 




^0" 


v 


e 


6 


7 


ep 




{kg.m"^) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 


(km.s""*^) 


(IcnLS^-*^) 


(GPa) 


(%) 


(%) 


(%) 


(%) 


HCP Iron 


Mao et al. fl998) 


12600 


1533 


367 


583 


1544 


835 


11.47 


5.92 




-0.36 


-F29.60 


-20.50 


2.46 


Laio et al. f2000l 


12885 


1697 


809 


444 


1799 


757 


11.45 


5.90 




-2.84 


-11.12 


+2.73 


0.80 


Vocadlo f 20071 


13155 


1311 


159 


1642 


1642 


1074 


11.14 


4.04 




+2M 


-15.22 


+15.72 


2.22 


BCC Iron 


Belonoshko et al. (20071 


13850 


1561 


1448 


365 






11.54 


4.20 


-617 








2.79 : 


Vocadlo f 20071 


13155 


1603 


1258 


256 






11.29 


4.11 


-167 








0.86 


Vocadlo f 20071 


13842 


1795 


1519 


323 






11.83 


4.24 


-370 








1.60 



Table 1 

Elastic properties of hexagonal (hep) and cubic (bcc) iron crystals and associated untextured aggregate: density (p), elastic constants 
(Cjj), hexagonal anisotropic parameters {e,5,j), cubic anisotropic parameter (z^), P-wave and ^-wave Voig t veloc ities (Vn^, Vj^'^), degree s 



of inhomogeneity {£p,^s)- The iron models have been computed at inner core temperature except 



Mao et al 



(|l998l l and 



Laio et al 



pood i 




Fig. 1. Modelling of the superficial layer at the top of the solid inner core as 
an untextured aggregate of iron "patches", with typical size d. ICB: Inner Core 
Boundary. 
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Fig. 2. P-wave velocity (a) and S'-wave velocities (b and c) in a single crystal of 
iron as a function o f prop agation direc tion with respect to symmetry axis, f or 3 iron 



models: 



Laio et al. 



(|200m (solid line), 



Mao et al. 



(11993) (dashed line) and 



Vocadlo 



(j2007l ) (dotted line). Q is the angle between the direction of propagation and the 



symmetry axis of the crystal. 



P waves 



S waves 



0,000 



-0,002 - 



-0,006 - 



-0,008 - 



-0,010 



— 1 1 — I I I I I I I 1 1 — I I I I I 1 1 1 1 — I I I I I I 




HOP 



J I I I I I I I 



0,00 



- -0,02 



- -0,04 



- -0,06 - 



- -0,08 - 



-0,10 



— 1 1 — I I I I I I 



— 1 — I I I I 1 1 1 1 1 — I I I I I r 



— Laio et al. (2000), T=OK 
Mao et al. (1999), T=300K 
Vocadio (2007), T=5500K 



HCP 



_i I 



J I I I 



10' 



10"' 10" 



10^ 10"^ 



10"' 10" 



10' 



— I — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 11 I 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 

0,000 1- 0,00- 



- -0,05 - 



-0,015 



-0,020 




BCC 



J I I I I I I I 



-0,10 



-0,15 



-0,20 



- — Belonoshko et al. (2007), T=6000 K. 
- ■ Vocadio (2007), T=5500K 
■ ■ ■ Vocadio (2007), T=6000K 



BCC 



J I I I I I I I 



10"' 



10"' 10" 

k^a 



10^ 10"^ 



10"' 10" 

k a 



10' 



Fig. 3. Normalized variations of longitudinal-wave (left) and transverse-wave (right) 
phase velocities in hexagonal iron aggregates (top) and in cubic iron agregates 
(bottom) as a function of adimensional frequency ko a. ko and Vq are the Voigt 
wave vector and velocity respectively, Ve is the effective velocity in the aggregate, 
and a is the correlation length. 




Fig. 4. P-wave attenuation (1/Qp) as a function of adimensional frequency koa. 
Amplitude attenuation is computed for three hexagonal iron models and three 
cubic iron models proposed for the inner core. The grey dashed lines correspond 
to Qp = 300 and Qp = 600. 
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Fig. 5. Number of acceptable iron crystals as a function of the correlation 
length: hexagonal symmetry (left) and cubic symmetry (right). The corresponding 
aggregates verify the seismic observations at 1 Hz: P-wave and S-wave velocities are 
in the range of the ak 135 model of the whole inner core, and attenuation quality 
factor Qp lies between 300 and 600. 




Fig. 6. Density of acceptable hexagonal crystals as a function of the anisotropic 
parameters e and 6, for a = 140 m (top), a = 220 m (middle) and 
a = 1000 m ( bottom). For compari s on, co lor dots r epresenting eigh t hey ir on models 



are plotte d: 



Mao et al 



(green) , 



Stixrude and Cohen 



199811 rbrown'l. 



I995II fredl. 



Steinle-Neumann et al 



Stein 



e-Neumann et al 



(I2OO1I) (purple). 



Soderlind et al 



19991'! (cyan 



Vocadlo et al 



1996ll (black) . 



;aio et al 



(I2OOOI) 



2003al ) (orange). 



Vocadld (|2007l ) (magenta). Only two models have been computed at inner core 



temperature (squares). 



a = 140 m Number of acceptable models = 9750 
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a = 220 m Number of acceptable models = 154345 
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a = 1000 m Number of acceptable models = 5199 
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Fig. 7. Distribution of Voigt velocities of acceptable hep iron crystals for different 
correlation lengths: a = 140 m, o = 220 m, a = 1000 m (from top to bottom). 
Vertical lines delimit the range of seismic velocities in the inner core. 



a = 150 m Number of acceptable models = 114 
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Fig. 8. Distribution of acceptable elastic constants Cij / p of cubic iron for different 
correlation lengths, indicated on top of each plot. At large correlation length, we 
observe two distinct families of cubic models (sho wn in gr e y and black). Circles 



correspond to two hcc iron models proposed by 



Vocadld (j2007l ) for a density 



p = 13155 kg.m ^ and a temperature T=5500 K (black), and for p = 13842 kg.m 
an d T=6000 K (whi t e). T 



by 



Belonoshko et al. 



re square corresponds to the hcc iron model proposed 
(|2007l ) for a density p = 13580 kg.m"'^ and a temperature 



T=6000 K. 
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Fig. 9. Distribution of Voigt velocities for the two families of acceptable cubic models 
of iron shown in Figure 8. The correlation length is indicated on top of each plot: 
a =150, 200, 500 m. Circles correspond to the two bcc i ron models propo s ed by 



(120071) 



Vocadld (j2007l ) and the square is the bcc model proposed bv lBelonoshko et al 
(same symbols as in Figure 8). Vertical lines delimit the range of seismic velocities 
in the inner core. 



